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■ We perform a QCD analysis for the radiative decay of a heavy ^Si quarko- 

nium into the tensor meson /2(1270). We make an attempt to separate the 
Q_i' nonperturbative effect related to the quarkonium and that related to the ten- 

iu: 



sor meson, the former is represented by NRQCD matrix elements, while the 
later is parameterized by distribution amplitudes of gluons in the tensor me- 



^ ■ son at twist-2 level and at twist-3 level. We find that at twist-2 level the 



helicity A of the tensor meson can be and 2 and the amplitude with A = 2 is 
suppressed. At twist-3 level the tensor meson can have A = 1. A comparison 
with experiment is made, an agreement of our results with experiment can 
be found. We also briefly discuss the radiative decay into tj and obtain a 
prediction for T — > 7 + r]. 

PACS numbers: 13.25.Gv, 14.40.-n, 12.38Bx. 
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1. Introduction 

Decays of heavy quarkonia into light hadrons are forbidden processes by OZI rule, their 
study will help us to understand how gluons, which are fundamentally dynamical freedoms of 
QCD, are converted into light hadrons. Radiative decays of quarkonia into a single hadron 
provide an ideal place to study the conversion of gluons into the light hadron. Because 
there is only one light hadron in the final state, there is no complication due to complicated 
interactions between light hadrons, and it is relatively easy to determine the conversion of 
gluons into the light hadron. In this work, we propose to study the radiative decay into 
/2(1270), which is a spin-2 particle and has the mass m = 1.27GeV. 

The radiative decay of J/\l/ was observed long time ago, and an analysis of the polariza- 
tion of /2 was also performed where the main decay mode /2 ^ n + n was used. Recently 
the decay of T was also observed by CLEO , the branching ratio was determined. 

In the decay the initial state is a heavy quarkonium, say T, which can be taken as a bound 
state of a b- and 6-quark as an approximation. Then the radiative decay of T undergoes 
as the following: the quarkonium will be annihilated into a real photon and gluons, the 
gluons will subsequently be converted into the tensor meson /2. In the heavy quark limit, 
the 6-and 5-quark moves with a small velocity v, hence an expansion in v can be employed, 
nonrelativistice QCD(NRQCD) can be used to describe the nonperturbative effect related to 
T P]. Also in the limit, the tensor meson /2 has a large momentum, this enables an expansion 
in twist to characterize the gluonic conversion into /2, the conversion is then described by 
a set of distribution amplitudes of gluons. The large momentum of /2 requires that the 
gluons should be hard, hence the emission of the gluons can be handled by perturbative 
theory. The above discussion implies that we may factorize the decay amplitude into three 
parts: the first part consists of matrix elements of NRQCD representing the nonperturbative 
effect related to T, the second part consists of some distribution amplitudes, which are 
for the gluonic conversion into /2, the third part consists of some coefficients, which can 
be calculated with perturbative theory for the 66-pair annihilated into gluons and a real 
photon. This work is an attempt to obtain a factorized form for the decay amplitude 
at tree level, in which we analyse the gluon conversion up to twist-3 level and give the 
definition of the corresponding distribution amplitudes. At loop levels, if the factorization 
still holds, the perturbative coefficients should be free from infrared singularities, and all 
nonperturbative effects should be contained in the NRQCD matrix elements and in the 
distribution amplitudes. To warranty this one needs to prove the factorization at loop 
levels, which is beyond the current work. 

The radiative decay was studied in the framework of a nonrelativistic quark model [Q, 
in which, not only a nonrelativistic description is employed for the initial quarkonium, but 
also for the tensor meson. The later is not well justified. Instead of this we use a set of 
gluonic distribution amplitudes for /2, which are relativistic, gauge invariant, and universal. 
Hence once information of them is extracted from one experiment, it can also be used for 
predictions for other experiments. For example, recently hard exclusive production of /2 was 
studied [Q, in which the same distribution amplitudes appeared. In this work we will use 
our results for the two decays: T ^ 7 + /2 and J/\l' ^ 7 + /2, and an agreement between 
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experiment and our results can be reached under certain conditions. 

Part of our results may also used for the radiative decay into rj. It turns out that the 
decay amplitude is at twist-4 level. Without a complete analysis we can still predict the 
decay width for T — 7 + We will briefly discuss this decay mode. 

Our work is organized as the following: In Sect. 2 we introduce our notation and perform 
the analysis at twist-2 level. In sect. 3 we perform the analysis at twist-3 level. In Sect. 
4 we make a comparison of our results with experiment and give a brief discussion for the 
decay into rj. Sect. 5 is our conclusion. 

Throughout of our work we take nonrelativistic normalization for the quarkonium state 
and for heavy quarks. 

2. Notations and results at order of twist-2 



We consider the exclusive decay of T in its rest-frame: 

T{P)^j{q) + Mk), (1) 

where the momenta are given in the brackets. Since /2 is a spin-2 particle, its polarization is 
described a symmetric, trace-less tensor e'^'^{k, A), where A is the helicity of f2- This tensor 

can be constructed by introducing polarization vectors cj's. We take a coordinate system in 
which /2 moves in the direction of the ^-axis. In this system the vectors take the form 

= ^(0,1,^,0), a;^(-l) = -L(o, 1,-^,0), a;^(0) = i-(|k|,0,0,A;°). (2) 

With these vectors and with Clebsch-Gorden coefficients one can construct e'^^{k,\) for 
A = 0, ±1 ± 2. At the leading order of QED the -S-matrix element for the decay is 

(7/2(A)|5|T) = -ieQ,e*^ ■ J dSe''^-^MX)\b{zh,b{^)\T), (3) 

where Qb is the electric charge of 6-quark in unit of e, b{x) is the Dirac field for 6-quark, e*^ 
is the polarization vector of the photon and A is the helicity of f2- In Eq.(3) we only take 
the part of 6-quark in the electric current into account. If two gluons are emitted by the b- 
or 6-quark, we obtain the corresponding contribution to the ^'-matrix element 

{h{X)\T [b{xh ■ G{x)b{x)b{y)^ ■ G{y)b{y)b{z)Yb{z)\ |T), (4) 

where G{x) is the gluon field. Using Wick-theorem we can calculate the T -ordered product 
and we only keep those terms in which one b-field and one &-ficld remain uncontracted. Then 
the matrix element takes a complicated form and can be written in a short notation: 

W2 = -%\eQ,9% I rfW|/dSd^a;iAie*''"(/2(A)|G'^(a;)G'^(|/)|0) 

{0\bj{x,My,)\T) ■ M^r'''\x,y,x,,y,,z), (5) 
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where Mj^'''''"^(x, 7/, xi,yi, z) is a known function, i and j stand for Dirac- and color indices, 
a and b is the color of gluon field. The above equation can be generalized to emission 
of arbitrary number n of gluons, the corresponding contribution is Wn, then the ^'-matrix 
element is the sum 

(7/2(A)|5|T)=^«;„. (6) 

n 

In each contribution there is the same matrix element {0\bj{x)bi{y)\T) . For this matrix 
element the expansion in v can be now performed, the result is: 

{0\b,{xMy)\T) = -^(P+yi^-).,(0|xVV|T)e-^-(-+^) + 0{v'), (7) 

where x^('0) is the NRQCD field for b(b) quark and 

P± = (l±7°)/2, 

p'^^K, 0,0,0), (8) 

where is the pole-mass of the 6-quark. The leading order of the matrix element is 
0{v^), we will neglect the contribution from higher orders and the momentum of T is then 
approximated by 2p. It should be noted that effects at higher order of v can be added with 
the expansion in Eq.(7). Taking the result in Eq.(7) we can write the ^'-matrix element as: 

W2 = t^eQ,gl{27r)'6\2p - k - g)£;(0|xVV|T) 

/ ^r2^-(^'?i'^)-^r'(p,A:,?i), (9) 

with 

rf (A;,gi,A) = / dx^e-«-^(/2(A)|G'«''^(a;)G"^'''(0)|0). (10) 



To obtain Eq.(9) we have used the color-symmetry and the translational covariance. 
^2'{k^Qi) contains all nonperturbative effect related to /2, while R2'^''^{p,k,qi) is a per- 
turbative function, whose physical interpretation is that it is the amplitude for a ^5*1 bb pair 
emitting two gluons and a real photon, and the quarks have the same momentum p. The 
contribution of W2 may be represented by the diagrams given in Fig. 1. In Fig. 2 we give the 
diagrams corresponding to emission of three gluons, their contributions will be analysed 
in the next section. Emissions of more than three gluons will lead to contributions which 
are at orders higher than those of twist-3, and will not be considered in this work. 

To perform the twist expansion it is convenient to introduce the light-cone coordinate 
system, in which a vector A is given by A'^ = (A"*", A^, A^, A"^) = (^4+, A^, At), the compo- 
nent A'^ and A~ is related to A^ and A^ in the usual coordinate system by 

A-^^^iA' + A% A-^^{A'-A% (11) 
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In the light-cone coordinate system the momentum k of /2 takes the form: 

2 

Ttl 

F = (fc+,A;-,0,0), k- = ^. (12) 

In the heavy quark limit, is very large, while k~ goes to zero. We also define two vectors 
n and / and a tensor dr- 

n^ = (0,l,0,0), r = (1,0,0,0), 

= gt^'^ - n^r - Tfr. (13) 

Different components of a vector can be projected out with these vectors and the tensor: 

A+ =n-A, A- = l-A, A'i. = di^^A^. (14) 
The light-cone gauge is chosen in our work and is defined by 

n.G(x) = G+(x) = 0. (15) 
In this gauge a set of components of the gluon field strength takes a simple form: 

G+^'ix) = d+G^'ix) = -^Cix). (16) 

The nonperturbative object Fg'^ characterize the conversion of two gluons into f2- The x- 
dependence of the matrix element {f2{k, A)|G'°'^(x)G"''^(0)|0) is controlled by different scales: 
the ^"-dependence is controlled by k~^, while the x~^- and xy-dependence is controlled by 
the scale Aqcd or k~ , which are small in comparison with k~^. Because of these small scales 
we can expand the matrix element in x+ and in xt- With this expansion we obtain 

rf (A;,gi,A) = {27rf 6 {q^)6\qrT) ■ J rfx-e'^^^" (/2(A) |G"''^(a;-)G"'^(0) |0) 

+t{27rf5{q^)^5''{q,T)- J rfx-e'^^^^" (/2(A) |9^G"''^(x-)G"'^(0) |0) 
+^{27rf--5iq^)5^{q^T)■ f rfx-e'^''^" (/2(A)|9-G'^''^(x-)G"^'^(0)|0) 



dqi 

+ , (17) 

where we introduced a short notation: 

G^'ix-) = Cix-n). (18) 

In Eq.(17), the leading twist of the first term is 2, while the leading twist of the second term 
is 3. By considering Lorentz covariance and varying the vector n |^ one can show that the 
derivative d~ in the third term is related to d^, hence the leading twist of the third term is 

4. The stand for contributions which have twist more than 3 and will be neglected. 

It should be noted that the expansion in Eq.(17) is equivalent to the coUinear expansion 
around = {qf, 0, 0, 0) for R'^^'^^ip, k, qi). 
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The leading-twist contribution in the first term is specified by that the indices jj, and u 
are all transversal indices, i.e., or z/ = 1,2. By considering the covariance under Lorentz 
boost along the ^-axis and the invariance of rotations around the 2;-axis we find at twist-2 
for A = 0: 

rr(^, q^, 0) = i2nr5{qT)5\q^T) ■ 2k+x,{x, - 1) "^^^ ^^^^ 
where qi — Xik~^, Fq{xi) is the gluonic distribution amplitude for A = and is given by 

i^o(xi) = ^J dx-e-^^^'^'^-{MO)\G^'^^{x-)G-\{Om (20) 
It should be noted that 

-^4^^ = eno) = elUm'T'd'^;^. (21) 

At twist-2 level r^^''{k, q-i, A) can also be nonzero for A = ±2. The contribution for A = 2 
can be written as 

rr (^, gi, 2) = {27:)%q^)5\q,T) ■ ^^^^^^^ _ ^^ g*^'^(2)F2(a:i), (22) 
where F2{xi) is the gluonic distribution amplitude for A = 2 and is defined by: 

F2{x{) = ^ J c^x-e-^^i'=^^"(/2(2)|G"'+'^(x-)G'»'+^(0)|0)a;^(l)u;,(l) (23) 

The polarization vector s^'^{2) is given by 

£'^'^(2) =u;'^(l)a;'^(l). (24) 

At twist-2 level /2 can not have the hehcity A = 1, because the two gluons are coUinear 
with a relative angular momentum which is zero and they have Bose-symmetry. With these 
results for ^^^{k, qi, A) and with the exphcit form of B^^^^lp, k, qi) we obtain the -S-matrix 
element with A = and with A = 2: 

(7/2(0)|5|T) = '-eQ,g',{2n)'S\2p - k - q)e*'{Q\x^a'^P\T) ■ ^ ■ To, 
D 

(7/2(2)|5|T) = '-eQ,gli2nr5\2p - k - g)(-£>*'^(l))u;*^(l)(0|xVV|T) -^-T,, (25) 



with 



/ 8ml{x^-l)+m\l-2x^) 
°" ' '2x,(l-x,)(Aml(x,-l)+m^(l-2x,)) 



/TYl 
"^^'xAl - x,)Umi(x, -1) + mHl - 2x,)) ' ^^^^ 



In the calculation we have kept the mass of /2. The effect of the mass is usually regarded 
as a correction to the leading twist effect. From the above results we can see that the 
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amplitude with A = 2 is suppressed by {m/mi))'^, although it is a leading twist contribution. 
At twist-2 two gluonic distribution amplitudes characterize the gluon conversion into /2 with 
the helicity A = and A = 2, respectively. They are defined in the light-cone gauge and 
are invariant under gauge transformations, which respect the gauge condition in Eq.(15). 
In other gauges a gauge link between the two gluonic operators in Fq and in F2 should 
be added to make them gauge invariant. Because of the momentum-conservation in the 
+- direction, the ftmctions Fq(x) and F2{x) becomes zero if ,t > 1 or x < 0. These functions 
are defined at certain renormalization scale fi, hence they depends not only on x but also on 
fi. The evolution of these twist-2 operators are well known. Our results given above receive 
corrections from orders of higher twist and the corrections are suppressed by {A/rrib)^ where 
A can be Aqcd or k~. 

With these results in Eq.(25) and Eq.(26) we complete the analysis at twits-2 level. 
The contributions at twist-3 come from the first and the second term in Eq.(17) and also 
from contributions with emission of three gluons, which are represented by Fig.2. We will 
analyse these contributions in the next section and show that the amplitude at twist-3 level 
is nonzero for A = 1. 

3. Twist-3 Contributions 

In this section we calculate the contribution at order of twist-3. It turns out that at 
this order the decay amplitude is nonzero only for A = 1. We will neglect in this section 
the effect of m. The effect may be included and including it will results in complicated 
expressions. We start with Eq.(17). There are two terms which are twist-3 contributions, 
one term is specified by that one of the indices fi and z/ is — in the first term, while another 
term is specified by the indices jj, and to be all transversal in the second term. Examining 
their Lorentz structure we find it is nonzero for A = ±1. Hence we obtain the leading- twist 
contribution for A = 1 as 

J dx-e-'«^''"(/2(l)|n'^G'"'-(x-)G'^'"(0)+n"G'^'^(x-)G'"'-(0)|0) 

where G^"" = d^'G^ For r^'^'^"(fc, qi) we can define a distribution amplitude 

Goix) = ^/ dx-e-^^'^^^-{MX)KT^-{x-)GrmO)u;,{l). (28) 

It should be noted that G^ in the light cone gauge is not an independent dynamical freedom. 
Its relation to the gauge field Gj. and quark fields can be obtained by solving equation of 
motion of QCD. One can obtain that can be expressed by Df^Gxn color-charge 
current of quarks, where is the covariant derivative in the adjoint representation of 

su{3y. 



J- 2 

■p(i)/ii/ 
i 2 



{k,qi, 1) 
{k,qi) 



{k,qi] 
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[D'^U^d'^Sab + gsfabcG"^'', fora,6 = l,---,8. (29) 

Because of Dj^Gtim the contribution is a twist-3 contribution. We will not try to give the 
result for but leave it with the compact form. Expressing V^^^'^ {k^qi) with Gq{x) we 
obtain the contribution to the ^'-matrix element (7/2(1) |5'|T)o: 



(7/2(l)|^|T)o = ^eg,5,2(27r)V(2p-fc-g)e*-a;*(l)r(0|xVV|T)- ^ 



./dx^Go(.). (30) 

The Lorentz structure of the integral in r^^^^(k, qi) can be decomposed as: 

C'^-p = J dx-e-"'^'''{f2il)\dl^G^^{x-)G^''{0)\0) 

- k^q^\k)^ I dx-e-^'^t^-{Miw^Gr'{x-)Gr'^mo) 

= 4fe+(^^^'_fe)+ {K^< + ^^'^^ - 3a;*''dr)Xo(xi) 
+ ^(a;*X" + - Suj*^di^')Ki{xi) 

+ ^{uj*'d^T + - 3a;*"cif )i^2(a;i)}, (31) 

where the functions are given by 

Ko{x) = ^ j da;-e-^^'=^^-(/2(l)|a^G^'+^(a;-)G^'+^(0)|0)u;p(l), 
K,{x) = ^/ cia;-e-^"-(/2(l)|9^G^'+''(x-)G^'%(0)|0)a;,(l), 
K2{x) = da;-e-^^'=^^-(/2(l)|a^G^'%(a;-)G^'+''(0)|0)a;,(l), (32) 

and they have the following property: 

Ko{x) = -Ko{l - x), K,{x) = -K2{1 - x). (33) 
With these results we obtain the contribution to the ^'-matrix element from r'2^^^{k, qi): 

(7/2(l)|^|T)i = ^eQ5^?,?(27r) V(2p - k - q)e* ■ a;*(l)r(0|xV>|T) 



(7/2(l)|5|T)2 = '-eQ,g',{2nrS'{2p - k - q)e* ■ u;*(l)r(0|xVV|T) 
1 1 



In the above results Go{x) is gauge invariant, hence the contribution (7/2(1) |S'|T)o also does. 
But the function Kq{x) and Ki{x) is not gauge invariant, the contributions in Eq.(34) also 
are not gauge invariant. The gauge invariance can be achieved by considering the emission 
of three gluons, which we will consider in the below. 

The contribution from type of diagrams given in Fig. 2 a can be written after the expansion 
in Eq.(7): 

W3a = '-eQ,g!{2n)'5\2p - k - g)£;(0|xVV|T) 

I (2j|4(2^)4 -'^3amM2M3(fe^ Qu 52, A) • R^'/^'"^''\p, k, qi, 52), (35) 

with 

r^™(fc,gi,g2, A) = J da;^dy^e-^''-^e-^''-Va(.c(/2(A)|G"'^H^)G''^^(l/)^^''^^(0)|0). (36) 

In Eq.(35) the factor i/A comes from the trace over color indices, the three gluons forms a 
C-even state, hence the contribution is proportional to fabc- The factor 1/3! is a statistical 
factor. The leading order of w^a is twist-3, its contribution is given by keeping the x~- 
and y^dependence in the matrix element only and by all gluon fields being transversal. By 
examining the Lorentz structure of this contribution it is nonzero only for A = ±1. The 
leading term of T^a for A = 1 can be written: 

r^'^'ik, Qi, 12, 1) = {2nmq^)^{q^)^\q^T)S'{q2T)C^r'^ (37) 
where the tensor Cg^'*^^^ is given by 

Cg^^l/^^W ^ i(t^*M2^W« ^^*M3^M1M2 _ 3u;*^ld^'^')Ho{Xi,X2) 

+ ^(a;*^i<^^^ +a;*^3<i''^ - 3u*^''dt^''"')Hl{xl,X2) 
+ ^{u;*'''df^'^' + u*^^d^^^' - 3u;*'''di^'^')H2{xi,X2). (38) 
The functions are defined by: 
^^0(^1,^2) = / dx-dy-e--i'=^----^'=^^-/„6,(/2(l)|G^'^(x-)G'^''^(y-)G'^(0)|0)u;^(l), 
Hi{x,,X2) = I dx-dy-e-'-''^---'-'''^y-Uc{f2il)\G'''''{x-)G'^^i^^^^ 
Hs{x,,X2) = / dx-dy-e--i^^----^'=^^-/„;.c(/2(l)|G'"''^(x-)G't(y-)G'5^'^(0)|0)a;^(l), (39) 

where the variable Xi and X2 are related to qf and q2 by 

qf — xik'^ and = X2k'^, (40) 
respectively. The functions are related to each other: 

Ho{xi,X2) ^ -Hi{x2,Xi), H2{xi,X2) ^ Ho{l - Xi - X2,Xi), (41) 
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and there are also some identities among them, one of them, which will be used later, is 
given by 

H2iXi,X2)=-H2ix2,Xi). (42) 

After a tedious calculation we obtain the leading-twist contribution of T^^ to the S'-matrix 
element 

(7/2(l)|^|T)3 = '-eQ,g',(2n)'5\2p - k - q)s* ■ a;*(l)f (0|xVV|T) 
igs (k+y f, , xi + 2x2-1 



f dxidx2 — -^—o ■Ho(xi,X2). (43) 

J X^X'? — X^ + Xo — Xo 



\f2ml (27r)2 J X\X2 -X\^x\-X2 

We find that this contribution can be combined with (7/2(1)! 51 T)i to form a gauge-invariant 
contribution. For doing this, we insert an identity into (7/2(1) |5'|T)i: 



we obtain 



dy-e-'^'^^^y = j dx2S{x2) = 1 (44) 



(7/2(l)|-5|T)i = (• • I dx,dx2-^ 



1 - 2xi 



{2nyJ ' '{l-xifxl 
j da:-dy-e-^^i'=^^"-^^^'=^^"(/2(l)|9^G'^+''(x-)G'^+^(0)|0)6^p(l), (45) 

where (• • ■) stands for some common factors. For (7/2(1) |S'|T)3 we make the variable change 
xi — > X2 and X2 — > xi and write some gluon fields into field strength tensors: 

(7/2(l)|5|T)3 = (• • ■)^, I dx,dx2 + 2x1-1 



(27r)2 / Xi{l - Xi- X2){xiX2 - X2 + x\- Xi) 

j dx-dy-e-'^'^^^--'^-^^y-hUf2{l)\Gr{y-)G'^^''{x^^^ (46) 

where (■ ■ ■) stands for the same common factors as in Eq.(45). We note that by setting 
X2 = Q the function in the first fine of Eq.(46) is the same in the first line of Eq.(45). 
Because of the identity in Eq. (44) we can combine the two terms together as 

(7/2(l)|5|T)i + (7/2(l)|5|T)3 

. . 1 f , , X2 + 2xi - 1 

— (■ ■ V7 To / dXidX2 7 -7 7, r 

[2-kY j Xi[\ — Xi — X2)[XiX2 — X2 + xf — Xi) 

J dx-dy-e-''''''^''~-'''^''^y~{f2{l)\d^G'''+''{x-)G'''\{0) 
-gsfabcGr(y')G''-'''(x-)G'^'\{0)\0)u;^(l). (47) 

We note that the two terms in the matrix element of the above equation can be written as 
a one term with the covariant derivative D!^{y) which is defined as 

[D''{y)lb = d'^^ab + 9sfabcG''''{y), for a, 6 = 1, • • • , 8. (48) 
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With this observation we can define our second distribution amphtude Gi{xi,X2): 

G,{x,,x,) = ^/ dx-dy-e-^^^'^^-^^^'''y'{Ml)m{y-)G^%x-)r^^^^ 

(49) 

This amphtude is gauge invariant in the hght-cone gauge. In other gauges two gauge hnks 
must be supphed to ensure gauge invariance. With it the sum can be written: 

<t/2(1)|S|T), + (7/2(l)|S|T)3 

Now we turn to the contribution represented by Fig. 2b, where the three-gluon vertex is 
involved. This contribution can be obtained from W2 in Eq.(9) by using perturbative theory 
with the three-gluon vertex for r2 once. We denote the term by using the perturbative 
theory for r2 as r23, which is given by: 

-iD''''{q,)d,G;{x)G'''P{x)G'''''{0) + q,,D/{q,)G'''''{x)G'^''{^^^^ 

+{lJ, ^ly, n, qi^k- qx)} (51) 

where D^^^^qi) is the gluon propagator in the light-cone gauge: 

^.^M.4(,.^_!^!£l±I^). (52) 

r \ n-q ) 

The leading twist of this term is 3, whose contribution comes from the second term in Eq.(51) 
with dfj — — . Again, here the tensor meson can only have A = ±1. We denote the 
corresponding contribution to the ^'-matrix element as (7/2(1) |S'|T)4, which is then given 

by 

(7/2(l)|^|T)4 = i^^eQ,glCl^fb\2v - k - q)e%Q\x^a'^\T) 

^,Rr\p, k, q,){i2nr5{q^)5'{q,T)gsD-^{q,) 

I dx-e-^''^^-fabc{f2{l)\G'''^''{x-)G\{x-)G'^'''{Q)\Q) 

+{lJ, ^ly, u^fj,, qi^k- qi)}. (53) 

At the first glance the above contribution may be divergent because of the gluon propagator 
is singular if we set qi — {qi ,0,0,0). However we find that the integrand is finite when 
Qi ~^ (?i^) 0)0,0), and it leads to a finite contribution which is the same if we replace the 
gluon propagator by its special propagator D^^[qi): 

D'-'^iqi) - D'-'^iqi) = (54) 

[n ■ qiY 
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A discussion of the special propagator can be found in [0, where it is shown that by using 
the special propagator of quarks twist-4 contributions to the hadron structure functions can 
be conveniently analysed. 

By taking care of the singularity the calculation is straightforward. We obtain 

(7/2(l)|5|T)4 = '-eQ,gl{27^Y5\2p - k - q)e* ■ a.*(l)r(0|xVV|T) 



^9s 



y/2ml J xf(l-xi) 

with Jo{xi) defined as 



dxi^j- ■ Jo{xi) (55) 



1 



27^A;^ 



To ensure the gauge invariance of the S'-matrix element, the sum (7/2(1) |S'|T)2 + 
(7/2(1) |'S'|T)4 should be gauge invariant. As the expressions for both terms stand, it does 
not look like that the sum can be written in a gauge-invariant form. A re-arrangement is 
needed. We first write Jo{xi) as an double distribution amplitude by using 

Giix-) = ^J dx2dy-e-^-^^^'^y--^-Gl{y-). (57) 
Then the contribution can be written as 
(7/2(l)|5|T)4 = (---)77^ / dxrdx2 



(27r)2 y (xi + 0:2)^(1 - Xi - X2) 

rfx-rfi/-e-^^^^-'^--"^'=''^-/,,,(/2(l)|G^'+^(x-)G;:(y-)G^'t(0)|0)a;^(l), (58) 

where (■ ■ ■) stands for the same common factors as before. We add to (7/2(1) |S'|T)4 a term 

^ dxidx2j ^ rj^ -H2{xi,X2), (59) 



(27r)2 / (xi + a;2)(l - xi - X2 

which is identically zero because of the property in Eq.(42). With this term the contribution 
can be written: 

(7/2(l)|5'|T)4 = (■■OtJt^ j dxidx2f2ixi,X2) 



rfx-rfy-e-^^'^'="'^--^^^'=^^"/.,,(/2(l)|G^'+^(x-)G;:(l/-)G^'t(0)|0)c.^(l), 
^ / 2 2(1-2x1-20:2) ] 

= i(x, + X2)^(l-Xi-X2) - Xi(x, + X2)(l-X,-X2)4 ^^'^ 

By setting X2 = the function f2{xi,0) is the same as that in the front of -^'2(3^1) in Eq- 
(34). With the same trick used to derive Eq.(47) we can write the sum as 



(7/2(l)|5|T)2 + (7/2(l)|5|T)4 

, . 1 



dXidX2f2{Xl,X2) 



c/x-d?/-e-"^'=''^'"-^^2^^^"(/2(l)|5^G"'+^(x^)G'''+'^(0) 
+9sfabcG^'^''ix-)Gliy-)G'^^+^iO)\0)Lu^il). (61) 
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Again the two terms in the matrix element can be combined together with the covariant 
derivative in Eq. (48). We can define the third distribution amphtude 6^2(2^1, 2:2): 



(62) 

The sum can be written: 

(7/2(l)|^|T)2+(7/2(l)|^|T)4 = (---)^ / dx^dX2f2(x„X2)G,(x„X2) (63) 

Adding different contributions together we obtain the S'-matrix element for A = 1 at twist-3 
level: 



i 



(7/2(l)|^|T) = -eQt9l{27ryS\2p -k- q)e* ■ a;*(l)r(0|xVV|T)^ • T, 
6 mt 



Ti = — [ / dxi — — ^ -Go{xi)+i I dxidx2{fi{xi,X2)Gi{xi,X2) 

+2f2{x,,X2)G2{x,,X2))], (64) 

where the function /i is given by 

r I X .X2 +2X1-1 

h[xi,X2)^^ T- (65) 

a;i(l - xi- X2)[xiX2 - X2 + xf- Xi) 

With the results given in Eq.(25), Ea.(26) and Eq.(64) we completed the analysis for 
the decay T — 7 + /2, in which the nonperturbative effect related to T and that related 
to /2 are separated, these effects are parameterized by NRQCD matrix elements and the 
distribution amplitudes Fq, Fi and Gi{i — 0, 1, 2), respectively. The distribution amplitudes 
characterize how gluons are converted into the tensor meson /2. Their definitions are given 
in this and last section. They are invariant under Lorentz boosts in the ^-direction, and are 
gauge invariant. With the analysis we find that the decay amplitude with A = 1 and A = 2 
are suppressed by the power A/m^ and (A/m^)^ respectively, in comparison with the decay 
amplitude with A = 0. 



4. Comparison with Experiment 



In this section we will compare our results with experiment and we will also discuss the 
radiative decay into rj. With the decay amplitudes derived in the last two sections we obtain 
the decay width: 

r(T ^ 7 + /2) = l^Qlaalim,) (l - " ^^C^lOl eS^)\T) Y,JT,\' (66) 

where we used 
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m^a'x\0){0\x^a^m = {T\Oj CS^)\T)e\e^y . (67) 

In the above equation e is the polarization vector of T, the matrix element (T|0j^(^5'i)|T) 
is defined in and the average over the spin is implied in the matrix element. As discussed 
in the last section the relative order of magnitude of the three amplitudes is given by 

|To| : \T,\ : {T^l = 0{1) : 0{—) : O(^), (68) 

where A = Aqcd or k~, and each amplitude receives correction at the order 0{^) relatively 
to its leading-order contribution. The corrections are unknown. For consistency we neglect 
terms with Ti and T2 and the mass m of /2 also should be neglected. Hence we approximate 
the decay width by 

r(T ^ 7 + /2) = In^Qlaa'Am,) (l - " :^^{mICS,)\r)\To\' 

.\l + Oi^) + 0{a,) + 0{v')], (69) 

where the orders of possible corrections are given in {■ ■ ■}. For T it is a good approximation 
to neglect these corrections, because rrib is large, while for J/\E' the corrections may be 
significant with the relatively small m^, each correction in {■ ■ ■} can be at order of 30%. We 
also neglect these for J/\E'. To compare with experiment we build the ratio 

Rh{T) = r(T ^ 7 + /s) 

r(T ^ light hadrons) 

27 2 " [ '^^ \ ^ \rp \2 \ 

where we used the result at leading orders for r(T ^ light hadrons), where light hadrons 
are produced through gluon emission. Because Feyman diagrams for r(T — »• light hadrons) 
have a similar structure as those given in Fig. 1, corrections from higher orders of and 
of f ^ can be cancelled at certain level in this ratio, we can expect that the ratio receives 
smaller corrections than the width does. This ratio does not depend on the nonperturbative 
property of T in our approximation. The value of Tq is unknown, it can be extracted from 
Rh{T), and used to predict Rh{J /^). As an estimation we neglect the //-dependence of Tq. 
With experimental data for hadronic branching ratios we obtain the ratio: 

^ fi(T^7 + /2) 



B{j/m ^ 7 + /2) 

0.923 Qlml a,(m,) ^ " 
0.707 ■ Qlml ■ a,(m,) ' 1 - 

0.059 (71) 



where we have used the parameters: 
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rUc = l.SGeV, mb = S.OGeV, a^lmj = 0.3, a,(m6) = 0.18. (72) 

Experimentally, the branching ratio has been measured recently, its value is [^: 

5(T ^ 7 + /2) = (8.1 ± 2.3+1?) X 10-^ (73) 

With the value for B[J /'^ ^ 7 + /2) = 1-38 x 10^'^ [|1^] we obtain the ratio determined by 
experiment: 

R ^ 0.058 ± 0.016+°:°^J. (74) 

This value is in agreement with our prediction. However, it should be kept in mind that 
experimental errors are large and corrections to our prediction can be substantial. Another 
uncertainty is from the pole-mass of c-quark, which value is not well known as the pole-mass 



of 6-quark. A recent QCD sum rule analysis for charmonium systems shows that rric 
can be 1.7GeV. If we vary the mass from 1.3GeV to 1.7GeV, then the ratio R in Eq.(71) 
changes from 0.044 to 0.076. 

As discussed before, we should neglect the decay amplitude with A = ±1 and with 
A = ±2 to consistently predict decay widths in our approximation. However the helicity 
amplitudes can be measured by measuring the polarization of /2, in which the main decay 
mode /2 TTTT can be used. This type of measurements have been done for J/\l/-decay 
|T[]. But, as pointed out in [^, the fitting formula used in the measurement for the angular 
distribution of 7 and tt was not correct. Hence we still do not know about various helicity 
amplitudes. With the sample of 5 x 10^ J/'ip^s, whose collection will be ended this year at 
BES, a new analysis will provide such information. In our approach, the angular distribution 
will be approximated by 

-j^, (X {(1 + cos2^,)(3cos2 0. - 1)W 

a cos 6^d(pT,d cos t^^ 

sin 26^ cos 6^, sin 20^(3 cos^ - l)Re(Tor*)} 
•{l + 0(^)} (75) 

where 6^ is the polar angle between 7 and the e+e~-beam axis, and 0^ are the polar and 
azimuthal angles of a pion in the /2-rest frame with respect to the 7-direction, and 0,r = 
is defined by the e"'"e~-beam axis. A general formula for the distribution is provided in [Q. 

With the experimental result in Eq.(73) one can estimate a phenomenological constant 
for the gluon conversion into f2- For this we write the distribution amplitude as 



i'o(x) = /f -/(a;), with / f{x) = l (76) 



where the constant fg has dimension 1 in mass and is defined by |^ 



(/2(fc)|G^^(0)G'^,(0)|0) = /J \ {{Kk, - \rr?go.,)e% - (« ^ /3)] - (/i 



{\9o.^e% - (« - m -(/"-- y)) ■ (77) 
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If we take the renormalization scale = as a very large scale, we may take the asymptotic 
form for the function f{x) at fi = rrib'- 



/(x) = 15x^(1 (78) 

with the experimental data we obtain the estimation ai fi = m^: 

^ 44MeV. (79) 

Before ending this section we briefly discussed the decay 

r^^{q)+vik) (80) 

This decay has been studied in [Q, based on a static quark model for rj, and the corresponding 

decay of J/ip also has been studied with various approaches P-[TT[]. In these approaches 



the decay is controlled by the f/A(l) anomaly, in which the conversion of gluons into 7] is 
characterized by a local operator of gluons. In our approach the corresponding operator 
becomes non-local. Following the results presented in Sect. 2 we can obtain the decay 
amplitude at twist-2 level: 



(7r/|5|T) = -leQ,g!{27ry6\2p - k - g)e,(0|xVV|T) e'^ ^ 

"'6 



ml f 1-2x1 



/ dx,—- ^F,(xi), (81) 



where the tensor e'^'^ is defined by 



e' 



t^^'-Xnp, (82) 



where e'^'^"^ is totally anti-symmetric with e^^"^^ = 1. The distribution amplitude F^{xi) 
characterizes the conversion of two gluons into 1], which is defined in the light-cone gauge as 

F,{x,) = ^ J (ix-e-^^i'=^^"(r7(A;)|G'^'+'^(x-)G"^'+^(0)|0)e^,. (83) 

Our result shows that the twist-2 contribution is suppressed by m^. In the twist-expansion 
the light hadron mass should be taken as a small scale as Aqcd, hence the contribution 
is proportional to Ag^^. This implies that a complete analysis should include not only this 
contribution but also twist-4 contributions, in which one needs to consider the contributions 
from emission of 2, 3 and 4 gluons. However, without a complete analysis we can always 
write the result of a complete analysis as 



(7r/|5|T) = eQ,g!{2nY6\2p - k - g)£,(0|xVV|T) e^"— ■ g„ (84) 

4o Ulfj 



where the parameter Qr, has a dimension 3 in mass. The parameter characterized the con- 
version of gluons into 7] and it does not depend on properties of T. Hence the above result 
also applies for the J/\l/ decay. Following the above analysis for /2 we obtain for 77: 
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g(T^7 + ^) 
5(J/^^7 + r/) 

0.923 Q^m^ a,(mj ^ - 



4m; 



0.707 Qlml as{mb) i - < 

^ 4.1 X 10"^ (85) 
With the experimental data for i?( J/\& ^ 7 + r^) we obtain 

5(T ^7 + ?]) ^ 3.5 X 10-^. (86) 
The predicted S(T ^ 7 + 77) may be too small to be observed. 

5. Conclusion 

In this work we have studied the radiative decay of a ^Si quarkonium into a tensor 
meson /2. A factorization at tree- level is performed for the decay amplitude, in which 
nonperturbative effect related to the quarkonium and that related to /2 is separated. These 
effects are parameterized with NRQCD matrix elements and with a set of gluonic distribution 
amplitudes, which are defined in this work. These amplitudes are gauge invariant and 
universal. At twist-2 level, we find two amplitudes characterizing the conversion of two 
gluons into /2 with the helicity A = and A = 2, respectively, and the contribution to the 
decay amplitude with A = 2 is suppressed by m^. At twist-3 level the decay amplitude is 
nonzero only for A = 1, three gluonic distribution amplitudes are introduced to describe 
the gluon conversion. At loops-level, the factorization may still be performed as shown in 
studies of other exclusive processes WM, but we will not study this subject in this work and 



leave it for future. 

In our approach an expansion in v, the velocity of the 6-quark inside T in its rest-frame, 
is used. We have only taken the leading order contributions at At this order, T can 
be considered as a bound state of the bb pair in a color singlet. The corrections to the 
leading-order results may also added. However, problems arise at order of v^. At this order 
T has a component in which the bb pair is in color octet and the component is a bound state 
of the bb pair with soft gluons. It is unclear how to add corrections from this component at 
order of v"^. It deserves a further study of these problems to understand a bound state of 
many dynamical freedoms of QCD. 

The NRQCD matrix elements are known from other experiment and form lattice QCD 
calculations, while there is not many information about the 5 gluonic distribution ampli- 
tudes. Using experiment data and our result we extract a parameter relevant for the A = 
amplitude. It should be emphasized that these distribution amplitudes are universal and can 
be used for production of /2 in other processes, like the one studied in ||^. In experiment, 
a measurement of the polarization of /2 can provide information about these distribution 
amplitudes with different A. With the sample of 5 x 10^ ^/^'s, whose collection will be 
ended this year at BES, such measurement is planed [|T^. 

With our results we can predict the ratio of the branching ratios for T and for J/\E'. An 
agreement between the prediction and experiment can be obtained. In the approach used 
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here, we can also perform an analysis for the radiative decay into r]. It turns out that twist-4 
effect should be included in a complete analysis. But without the complete analysis we are 
still able to predict the branching ratio for T and it may be too small to be observed. 
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Fig.l 

Fig.l: One of the diagrams for contributions of two-gluon emission. 
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(a) 




(b) 
Fig.2 

Fig. 2: Typical diagrams for contributions with three-gluon emission. 
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